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The crystal structure of a disordered form of Cri/3NbS2 has been characterized using diffraction 
and inelastic scattering of synchrotron radiation. In contrast to the previously reported symmetry 
(P6322), the crystal can be described by a regular twinning of an average P63 structure with three 
disordered positions of the Cr ions. Short-range correlations of the occupational disorder result in 
a quite intense and structured diffuse scattering; a static nature of the disorder was unambiguously 
attributed by the inelastic x-ray scattering. The diffuse scattering has been modeled using a reverse 
Monte-Carlo algorithm assuming a disorder of the Cr sub-lattice only. The observed correlated 
disorder of the Cr sub-lattice reduces the temperature of the magnetic ordering from 130 K to 
88 K and drastically modifies the field dependence of the magnetization as it is evidenced by the 
SQUID magnetometery. We conclude, that in contrast to the helicoidal spin structure assumed for 
P6322 form, the compound under study is ferromagneticafly ordered with a pronounced in-plane 
anisotropy. 

PACS numbers: 75.40.-s, 75.25.+Z, 61.05.cp, 61.05.fg, 75.50.Bb 


I. INTRODUCTION 

A lack of inversion symmetry in magnetic materi¬ 
als makes possible the antisymmetric Dzyaloshinsky- 
Moriya interaction (DMI) and results in a variety of mag¬ 
netic phenomena. A chiral helicoidal orderl^^, skyrmion 
lattice^^, a quantum blue foj^, a topological Hall 
effedP, a first order phase transition driven by critical 
fluctuation^ exemplify some of the topics attracting the 
attention of a wide range of physicists. An interesting 
and quite unusual circumstance is that, until recently, the 
most of the above listed phenomena have been proposed 
on the basis of the properties of MnSi and its substituted 
analogues. These compounds have a simple cubic crystal 
structure (P2i3 space group) with a unit cell which con¬ 
tains four metal and four silicon atoms in the 4a Wyckoff 
positions. The crystal structure is chiral in a sense that 
it does not coincide with its mirror image; for the ground 
state, and a certain part of the {H — T) (magnetic field 
- temperature) phase diagram, the magnetic structure is 
also chiral and its chi rality is unambiguously bound with 
the structural on#*2l. Recently, the class of chiral mag¬ 
nets has been augmented with monogermanides of some 
3d metal^i^*^ and Cu20Se0ji2®31, both having the same 
P2i3 space group as MnSi. 

Another interesting example is the compound 
Cri/ 3 NbSiP^ where a chiral soliton lattice has been pro¬ 
posed for its magnetic structur^i^. Its hexagonal struc¬ 
ture is built up from NbS 2 layers intercalated by Cr ions. 
There are three basic magnetic interactions which con¬ 
stitute the spin structure: (i) the ferromagnetic within 
the Cr layers Jl, (ii) another ferromagnetic one Jy and 


(iii) the DMI between Cr ions, the latter two interac¬ 
tions belong to the two intercalating layers separated by 
NbS 2 l^. The competition between the latter two interac¬ 
tions forms a helicoidal structure. Monte-Carlo simula¬ 
tions have recently showrP^ that the Cr ions are strongly 
ferromagnetically coupled in the (a6)-plane with J± ~ 74 
K, while they are weakly correlated between the layers 
with J|| ~ 10 K and D ^ 1.5 K. The first constant Jl de¬ 
termines solely the critical temperature of the compound, 
while the other two constants D and Jy balance in the 
length of the helix pitch. 

At variance with the crystallographically well studied 
B20 silicideJi^, the structural information on the other 
chiral magnets is rather incomplete and often limited to 
powder diffraction data. Thus, the absolute (chiral) sense 
has only been reported recently for both magnetism and 
crystal structure for Cu 2 0 Se 02 n^, this information is not 
yet available for Cri/ 3 NbS 2 because the chirality can¬ 
not be characterized from powder diffraction. As far as 
we are aware, the atomic coordinates in Cri/ 3 NbS 2 were 
measured only once, the refinement was done using nine 
neutron powder peaks and it was reported without any 
standard deviation^^. Other experimental structural in¬ 
formation is limited to the unit cell dimensionJ^^, and 
it is augmented with atomic coordinates borrowed from 
similar compound^^. In spite of the great interest in 
unusual and potentially technologically important prop¬ 
erties, the understanding of a microscopic nature of the 
chiral magnetism calls for more complete and precise in¬ 
formation on the underlaying crystal structure. 

The compounds Ma;NbS 2 (where M is a transition 
metal) are particularly interesting due to their layered 
hexagonal structure intercalated with transition metal 
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ions. Layered NbS 2 -based intercalates have revealed 
many fascinating phenomena including superconductiv¬ 
ity together with charge a nd sp in density wave^SMH^ 
Fermi surface nesting effect ^23111^ and also a structural 
disorder due to a distribution of the intercalated ions over 
several positions. Ma;NbS 2 materials also show a variety 
of magnetic properties as a function of intercalated M 
and x\ 

• paramagnetic for intercalated 

• previously reported as paramagnetic for 

later on as ferromagnetic below Tc = 50 K and 
another unidentified magnetic phase below T = 20 
K with X = 1/;P^ 

• antiferromagnetic for Co with Tjv = 26 K and for 
Ni with Tat = 90 K with a; = 1/pH^ 

• for Fe antiferromagnetic below Tm = 150 K with 
x = 1/4 an d a sp in-glass state below Tg = 40 K 
with X = I/ 322 IM 

• ferromagnetic and possible helicoidal state for Mn 
with X = 1/3 below Tc = 42 and with x = 1/4 
below Tc = 80 kPM 

• helicoidal for Cr with x = 1/3 below Tc = 130 
j dTsl^ 

In spite of the interest in the magnetic properties in 
Cri/aNbSiP^, there is an unexpectedly large spread in 
the reported temperatures of the magnetic ordering; the 
published values vary from lOCP^ to 132 iP^. 

Here we revisit the intercalated crystal structure of 
the chiral helimagnet Cri/ 3 NbS 2 using synchrotron radi¬ 
ation. All of the crystals we tested have a crystal struc¬ 
ture which is indeed chiral but it is different from the 
one reported earlier. The crystal structure of Cri/ 3 NbS 2 
reveals a disorder in the Cr sub-lattice; the presence of 
disorder has also been confirmed in a separate diffuse 
scattering measurement. The nature of the observed dif¬ 
fuse component relates to short-range correlations and, 
as evidenced by the inelastic x-ray scattering experiment 
(IXS), shows the elastic peak to be the major contributor 
of the characteristic diffuse scattering. 

The magnetic susceptibility measured with SQUID 
magnetometery shows Tc « 88 K, such a low value points 
to a link between the structural disorder and magnetic 
properties. Magnetization as a function of both field 
and temperature also differs from the data presented 
earliei!^, and expected from the theoretical consideration 
of helicoidal magnetJ^^. Taken together, the new data 
call for a revision of theoretical models for chiral magnetic 
interactions to include effects of disorder and short range 
correlations in spatial distribution of magnetic centers. 
Moreover, the attribution of the given magnetic tran¬ 
sition to the stoichiometry and disorder pattern might 
imply the revision of the data for other intercalates. 


II. EXPERIMENTAL 

A. Sample preparation 

A polycrystalline Cr 3 / 3 NbS 2 sample was prepared by 
sintering of the mixture of starting components at 800° 
in vacuum. The Cri/ 3 NbS 2 single crystals were grown 
using chemical transport method in iodine atmosphere 
at the temperature gradient T = 950 — 800 °C. The 
obtained crystals had a natural faceting and a metal¬ 
lic luster. The chemical composition has been indepen¬ 
dently checked with the Energy Dispersive X-Ray Analy¬ 
sis (EDX), with two randomly selected crystals and prob¬ 
ing different spots. The results are very homogeneous 
with Cr/Nb ratio ^ 0.3. 

B. X-ray Scattering Experiments 

Single crystal Bragg diffraction data were collected 
at the room temperature using the PILATUS@SNBL 
diffractometer at the BMOIA end station of the Swiss- 
Norwegian Beamlines at the ESRF (Grenoble, France); 
the wavelength of the synchrotron radiation was set to 
0.68290 A. The data were collected with a single (/-scan 
with angular step of 0.1° in a shutter-free mode with the 
PILATUS2M detector. 

The raw data were processed with the SNBL Toolbox, 
the integral intensities were extracted from the frames 
with the CrysAlisPro softwar^^Il^ the crystal structure 
was solved with SHELXS and refined with SHELXlP^. 

Single crystal data on diffuse scattering have been mea¬ 
sured at room temperature using the PILATUS6M de¬ 
tector at the ID23 (ESRF, Grenoble, France) with 0.689 
A wavelength and 0.1° slicing. The maps of reciprocal 
space were recovered from the raw frames using in-house 
software with the orientation matrix as refined in GrysAl- 
isPro. The diffuse scattering data were modeled with a 
reverse Monte-Carlo algorithm as implemented in the in- 
house software using structural models found from the 
Bragg scattering experiment. 

The IXS measurement was carried out at the beamline 
ID28 at the ESRF. The spectrometer was operated at 
17.794 keV incident energy, providing an energy resolu¬ 
tion of 3.0 meV full-width-half-maximum (Si(999) reflec¬ 
tion). IXS scans were performed in transmission geome¬ 
try along selected directions in reciprocal space. Further 
details of the experimental set-up and the data treatment 
can be found elsewher^^. 


C. SQUID magnetometery 

The magnetization of a single crystal sample was mea¬ 
sured using SQUID magnetometer (Quantum Design 
MPMS-5S). Temperature dependent measurements were 
performed in the range of 2 — 300 K. The field depen¬ 
dent magnetization was obtained in a field range from 
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H = —5to5TatT = 6K. The magnetic field was 
applied in two different orientations relative to the crys¬ 
tallographic directions parallel to the c-axis and within 
the a6-plane perpendicular to c. The magnetic moment 
per Cr atom has been calculated from the mass and the 
stoichiometry of the crystal based on the EDX data. 


III. RESULTS 
A. Magnetic properties 

The magnetization curves at T = 6 K for i? || c and 
H J- c show an anisotropic behavior (Fig. which differs 
significantly from the one of the canonical CrNbSj^. For 
H \\ c the magnetization curve consists of two regimes: 
at low fields up to Lf « 0.1 T the magnetization in¬ 
creases rapidly from 0 to 0.5 fiB/Cr. At higher fields the 
magnetization grows smoothly up to 1.2^^ at 5 T. Fast 
saturation of magnetization for H within the a6-plane 
and slow increase for H \\ c are typical features of mag¬ 
netic materials with uniaxial easy-plane anisotropy. For 
H A- c the saturation field equals to 50 mT which is be¬ 
low the range reported previously (from 1.5 to 2 

kO#^. The saturation magnetic moment Ms = 1.3^b 
per Cr ion is considerably lower than the spin-only value 
of 3/iB per Cr^+ ions with S' = 3/2. It is most likely 
that the reduction of the saturation moment is caused 
by delocalization of the d-electrons of Cr^+. 

The coercivity is smaller than 1 mT for both orienta¬ 
tions. This value, however, is comparable with a residual 
magnetic field of the SQUID magnetometer and, there¬ 
fore, the coercivity is below the experimental error. The 
negligibly small remanent magnetization and coercivity 
for T c can be explained by a continous spin rotation 
in the basal plane towards the external field direction. 
Along the c-axis, the field induces a canting of the spins 
competing against the easy-plane anisotropy. 

For H \\ c a linear M{H) dependence would be ex¬ 
pected in uniaxial easy-plane ferromagnets. However, a 
non-linear M{H) dependence is observed for the studied 
crystals which suggests that not only uniaxial anisotropy 
has to be taken into account in order to describe the mag¬ 
netization behavior. Most likely, also the non-uniform 
distribution of the intercalated chromium ions plays an 
important role. 

The temperature dependent magnetization at small 
magnetic fields (2.5 to 100 mT) well below saturation 
shows for both field orientations magnetic order below 
approx. 90 K (Fig. [^. The increasing field leads to 
an increase of the magnetization, while the moment is 
smaller along the c-axis than within the a6-plane. This 
is also reflected by the anisotropic magnetization curves 
in Fig. [2 

Arrott plots have been used to eliminate the influence 
of the external magnetic field and to accurately deter¬ 
mine the ordering temperature. According to the mean 
field theorji^, a plot of versus H/M for different 
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FIG. 1 . Magnetization curves for H || ab and R || c at T = 
6 K. 


(a) 




FIG. 2 . Temperature dependent magnetization for (a) H || ab 
and (b) H \\ c for different magnetic fields. 
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FIG. 3. Arrot plots for H |j a&-plane according to (a) a mean- 
field model {P = 0.5, 7 = 1 ) and (b) a 3D Heisenberg model 
iP = 0.365, 7 = 1.336). 


temperatures should show a linear dependence in the 
large field range. However, all curves are non-linear, 
but show a downward curvature (Fig. Therefore, 

we modified the Arrott plot taking the critical expo¬ 
nents into accountPSI. A plot of versus 

is shown in Fig. with the critical exponents along a 
3D-Heisenberg model (/3 = 0.365, 7 = 1.336)|2il repre¬ 
senting all linear behavior at high fields. From extrapo¬ 
lation of the slopes we obtain an ordering temperature of 
Tg = 88 K which is significantly smaller than reported 
ir |i 5 | 32 | critical exponents fulfill the Widom scal¬ 

ing relation (5=1-1- 7 /,^^, which gives 5 = 4.697 for 
the experimental data when the field dependence of the 
magnetization is considered as M (x at T = Tg. 

This result agrees well with S = 4.797 according to the 
3D Heisenberg model in contrast to 5 = 3 which would 
hold for a mean-field model. 


TABLE I. Three possible structural models for Cri/ 3 NbS 2 . 
SG is the space group, WP is the Wykoff positions. 


Model 

SG 

Ri Rw 

Gr WP 

Gr/Nb Ratio 

I 

P 6322 

0.148 0.29 

2 c 

0.333, fixed by symmetry 

II 

P 6322 

0.044 0.14 

2 c, 26 

0.4 

III 

P 63 

0.015 0.044 

26, 2 a 

0.314 





Model I 


Model II 


Model III 


FIG. 4. Fragment of crystal structures for the models I, II, III. 
The occupancy of disordered Cr sites is encoded as a colored 
sector. 


B. Bragg diffraction and crystal structure 


Bragg data were indexed with a hexagonal cell (a = 
5.738(5) A, 6 = 5.738(5) A, c = 12.018(5) A) in a close 
agreement with the unit cell dimensions reported pre¬ 
viously. Three models have been tested for the crys¬ 
tal structure, as compared in Table |l] and Figure 
The model I was proposed for the crystal structure of 
Cri/ 3 NbS 2 with space group P 6322 and the atomic po¬ 
sitions set close to those reported in Ref. [T71 Cr ion 
occupies the 2c Wyckoff position. The refinement gives 
rather high i?-factors, and analysis of the residual elec¬ 
tron density reveals additional maxima not accounted by 
the model. 

Model H includes an additional maximum at the 2b 
Wyckoff position which indicate presence of a disordered 
Cr ion. Disordered M 3 / 3 NbS 2 (where M is Mn, Fe, Co, 
Ni, V) has been reported, but the additional maxima 
were located at 2a and 2d and the positions are found to 
be only slightly populatecP^. Refinement with the model 
H has an acceptable quality, with, however, a few warn- 
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TABLE II. Crystal data for P 63 (model III). 


Empirical formula 

Temperature 

Wavelength 

Crystal system, space group 
Unit cell dimensions 


Theta range for data collection 


Cro.swNbSa 
293 K 
0.68290 A 
hexagonal, P 63 
a = 5.738(5) A 
b = 5.738(5) A 
c = 12.018(5) A 
1.63° to 25.25° 


Limiting indices 


-7 <= h <= 7, 
—6 <= k <= 6 , 


-15 <= I <= 15 


Reflections collected / unique 
Completeness to theta = 25.25 


2050 / 466 [Rint 
98.8 


0.0198] 


ing signs: (i) analysis of a variance of reflections included 
into the refinement indicates a possible twinning; (ii) the 
refined value of Cr content is higher than the measured 
with EDX (the Cr/Nb ratio is expected to be 1/3, mea¬ 
sured is 0.3, refined is 0.4). 

Therefore, the Model III, which has the P 63 symme¬ 
try and the twinning law 010 100 OOl, has been pro¬ 
posed. The Cr ions are found to be disordered over three 
independent crystallographic positions. The refinement 
gives the lowest i?-factors and the refined composition is 
close to the expected and measured with the Cr/Nb ratio 
around 0.314. Corresponding crystal data and results of 
the refinement at the room temperature are summarized 
in Tables [I| and [III and structural parameters are available 
as crystallographic information file (CIF) in supplemen¬ 
tary materials. 



hOI 


hhl 


FIG. 5. Maps of the diffuse scattering for different layers in 
reciprocal space. 


C. Diffuse scattering and structural disorder 


Diffuse scattering has a form of a honeycomb struc¬ 
ture extended along c*, the walls of honeycomb cells are 
weakly modulated along this direction (Fig. i- Bragg 
nodes stay in the centers of honeycomb cells. Along 
c* the nodes are connected by the diffuse streaks which 
points to the presence of 2D structural defects parallel to 
the NbS 2 planes; they can be interpreted as the bound¬ 
aries between the twins. 



FIG. 6 . (a) The /ihl-section of reciprocal space; (b) {E, Q) 
intensity map for the region of interest indicated in (a). 


D. Inelastic x-ray scattering 

Energy transfer scans along a wall of the honeycomb is 
shown at Fig. The intensity of the elastic peak follows 
the intensity of the diffuse component, while no substan¬ 
tial inelastic signal can be seen beyond. That confirms es¬ 
sentially elastic character of the observed scattering and 
therefore static nature of the underlaying disorder. 


E. Monte-Carlo modeling of the disordered 
structure 

Diffuse scattering has been modeled using the Reverse 
Monte-Carlo method (RMC) based on the Model II and 
using only one partially occupied Cr site (application of 
the twinned Model III was difficult because the number 
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FIG. 7. Diffuse scattering patterns: (exp) symmetrized ex¬ 
perimental, (calc) resilting from the RMC, (diff) difference 
map. The same intensity scale is used. 

10% 25% 50% 75% 90% 


FIG. 8. Randomly chosen xy-sections of the model clusters 
with the given Gr content. 


of correlating parameters was too large). 

The RMC simulation was performed on a cluster 32 x 
32 X 32 atoms to profit the efficiency of the fast Fourier 
transform. In an initial configuration the lattice sites 
were randomly occupied with a predefined concentration. 
Each simulation step included: 

1 . a swap between randomly chosen vacancy and 
atom; 

2 . calculation of the scattering intensities together 
with adjustment of the background and scaling fac¬ 
tor (c* dependence was removed by the averaging); 

3. calculation of and acceptance or rejection of the 
swap. 

For the given cluster size, a reasonable convergence 
(Fig. § is achieved in ~ 10® steps. As c* correla¬ 
tions were neglected, inspection of any XY section of 
converged cluster is illustrative enough. Qualitatively, 
one could summarize the message as follows: vacan¬ 
cies/atoms are distributed in a way to minimize contacts 
between the neighbors of the same type, locally resulting 
in a loss of hexagonal packing (Fig. |^. 


IV. DISCUSSION AND CONCLUSION 

Our results, from both Bragg and diffuse scattering 
unambiguously indicate that we are dealing with a dis¬ 
ordered form of the title compound. Notably, the crystal 


structure of this chiral helimagnet has never been ac¬ 
curately probed before with a single crystal diffraction 
experiment; we cannot, therefore, surely state that the 
disordered structure we observe is a generic one. How¬ 
ever, the magnetic ordering temperature is reported in¬ 
consistently in the literature, which strongly suggests the 
necessity of a careful structural characterization in order 
to validate the real structure of the material under study. 

We show that the correct symmetry is indeed chiral, 
but it is lower than the reported before. The P 63 struc¬ 
ture is twinned and disordered with respect to the Cr sub¬ 
lattice, which provides most of resonant contribution. As 
a result, the Flack parameter of the entire structure is not 
well defined. In contrast to MnSPand Cu20Se0^Sl the 
absolute structure cannot be reliably determined from 
the crystals tested in the present study. 

The crystal structure of the new form is almost identi¬ 
cal to the previously reported ones in terms of the geom¬ 
etry of the NbS 2 fragment. The main difference to the 
previously proposed structural model (P 6322 ) is a dis¬ 
order in Cr layers. In contrast to the initially expected 
P 6322 structure that assumes only one position for the 
Cr ions, the experimentally observed P 63 symmetry of¬ 
fers three Wyckoff positions with independent occupan¬ 
cies. Unconstrained refinement of the occupancies results 
in Cro. 3 i 4 NbS 2 , that is reasonably close to the stoichio¬ 
metric composition. 

The disorder found in the magnetic Cr sub-lattice is 
not random but correlated and results in the honeycomb¬ 
like diffuse scattering shown to be essentially elastic by 
IXS. This pattern suggests a tendency towards the clus¬ 
tering of Cr ions in hexagonal fragments within the lay¬ 
ers. Interlayer correlations are manifested in the form of 
modulated diffuse rods propagating along c*. Such a spe¬ 
cific correlated disorder should strongly affect magnetic 
ordering, in particular, the ordering temperature; more 
experiments on diffuse scattering combined with magne¬ 
tization measurements are needed to clarify an interplay 
between the correlated disorder and magnetic properties. 

We conclude with a note on the necessity of further re¬ 
vision of crystal structure of the chiral magnets. A chiral 
long-range magnetic spiral (hundreds or even thousands 
of Angstroms) is a giant object and it should not neces¬ 
sarily recognize the chirality of the few Angstroms size 
crystal unit cell. This is the case in rare-earth magnets 
like Ho, where the hexagonal structure is achiral and cen- 
trosymmetric but the magnetic ordering is chiral and its 
chirality can be changed by an external impadP^. How¬ 
ever, in the case of chiral magnets with a chiral crys¬ 
tal structure, as it was shown for B 20 magnetj' '*®B 6 ( 
atomic arrangement in a unit cell unambiguously dic¬ 
tates how the giant magnetic spiral must twist. This chi¬ 
ral coupling can be mapped phenomenologically by the 
Dzyaloshinsky-Moriya interaction, a microscopic mecha¬ 
nism behind it is still to be uncovered. Thus, the same 
link has been seen for both metal^ and insulatori^^ indi¬ 
cating that it is not related to the conduction electrons. 
In the 3d-metal monogermanides the link changes its sign 
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as a function of che mical composition, that points out the 
role of Sd-electron^^^Ell and, possibly, geometric factors 
related to the ionic radiP^. The results we present here 
indicate that a disorder could also affect the chiral mag¬ 
netism and illustrate that not only the average structure 
but also local deviations from it may influence the mag¬ 
netic properties. 
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